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Abstract 

Two-component systems are conceived as macrogases, and the re- 
lated equation of state is expressed using the virial theorem for sub- 
systems, under the restriction of homeoidally striated density profiles. 
Explicit calculations are performed for a useful reference case and a 
few cases of astrophysical interest, both with and without truncation 
radius. Shallower density profiles are found to yield an equation of 
state, <j) = <j){y,m), characterized (for assigned values of the fractional 
mass, m = Mj/Mi) by the occurrence of two extremum points, a 
minimum and a maximum, as found in an earlier attempt. Steeper 
density profiles produce a similar equation of state, which implies that 
a special value of m is related to a critical curve where the above men- 
tioned extremum points reduce to a single horizontal inflexion point, 
and curves below the critical one show no extremum points. The sim- 
ilarity of the isofractional mass curves to van der Waals' isothermal 
curves, suggests the possibility of a phase transition in a bell-shaped 
region of the {Oycj)) plane, where the fractional truncation radius along 
a selected direction is y = Rj/Ri, and the fractional virial potential en- 
ergy is (p = (£'jj)vir/(-E'jj)vir- Further investigation is devoted to mass 
distributions described by Hernquist (1990) density profiles, for which 
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an additional relation can be used to represent a sample of A'^ = 16 
elliptical galaxies (EGs) on the (Oycf)) plane. Even if the evolution of 
elliptical galaxies and their hosting dark matter (DM) haloes, in the 
light of the model, has been characterized by equal fractional mass, 
m, and equal scaled truncation radius, or concentration, = Ru/rj^, 
u = still it cannot be considered as strictly homologous, due to 
different values of fractional truncation radii, y, or fractional scaling 
radii, = rj/rj, deduced from sample objects. 
keywords - galaxies: evolution - dark matter: haloes. 

1 Introduction 

Ordinary fluids (e.g., gases and liquids) may be bounded by rigid walls which 
allow particle number conservation, avoiding evaporation. Macroscopical 
parameters (pressure, density, and temperature) remain uniform within the 
box, due to its reduced dimensions. On the other hand, astrophysical fluids 
(e.g., stars and galaxies) may be conceived as bounded by "gravitational" 
walls which violate particle number conservation by allowing evaporation. 
The macroscopical parameters exhibit gradients because they cannot remain 
uniform within the "gravitational" box, due to its large-scale dimensions. 

In particular, sufficiently extended celestial objects show at least two dis- 
tinct components: core-envelope for stars, core-halo for elliptical galaxies, 
bulge-disk for spiral and lenticular galaxies, baryonic-nonbaryonic for viri- 
alized (matter) density perturbations, and matter-dark energy for virialized 
(matter -|- dark energy) density perturbations. On this basis, an investigation 
on two-component astrophysical fluids appears useful for the comprehension 
and the interpretation of what is inferred from observations. To this aim, 
the choice of the density profiles is of basic importance. The laws of ideal 
and real gases were deduced from ordinary fluids, characterized by uniform 
density proflles. Accordingly, it is expected that astrophysical fluid laws are 
related to the specifled density profiles, and different laws hold for different 
matter distributions. 

Strictly speaking, a density profile should be deduced from the distribu- 
tion function, or vice versa. Unfortunately, the determination of the distribu- 
tion function is much more difficult than in one-component systems, and only 
a few cases have been studied in detail at present (e.g., Ciotti, 1996, 1999). 
On the other hand, global properties exhibited by simple density profiles 
(with somewhere negative distribution function) are expected to maintain 
a similar trend in dealing with much more complex density profiles (with 
nonnegative distribution function). 
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As the current attempt is mainly aimed to explore global properties in- 
stead of local properties, density profiles shall be selected according to their 
intrinsic simplicity, regardless from the physical meaning of the distribution 
function. Configurations described by simple density profiles could be suffi- 
ciently close to their counterparts described by self-consistent density profiles, 
and related results hold to a first extent. In any case, the self-consistency 
of density profiles with respect to nonnegativity of the distribution function, 
can be checked using a specific theorem (Ciotti and Pellegrini, 1992). 

To a first extent, the particle number shall be supposed to be conserved, 
which is equivalent to conceive the boundary of each subsystem as a perfectly 
reflecting surface in order to avoid evaporation. In this view, a possible choice 
of macroscopical parameters is: the fractional virial potential energy, 0; the 
fractional truncation radius, y\ and the fractional mass, m; as done in a 
pioneering paper with respect to uniform density profiles i.e. homogeneous 
configurations (Caimmi and Secco, 1990, hereafter quoted as CS90). Ac- 
cordingly, each subsystem is supposed to be virialized, in the sense that the 
virial equations are satisfied by averaging over a sufficiently long time, and 
particles move within a bounded region (e.g.. Landau and Lifchitz, 1966, 
Chap. II, §10; Caimmi, 2007a). Then virial equilibrium is a necessary (but 
not sufficient) condition for dynamical or hydrostatic equilibrium, which, on 
the other hand, does not imply pressureless configurations, as the stress ten- 
sor is related to the kinetic-energy tensor (e.g., Binney and Tremaine, 1987, 
Chap. 4, §2). 

For sake of simplicity, the applications of the general theory shall be re- 
stricted to homeoidally striated density profiles (e.g., Roberts, 1962; Caimmi, 
1993; Caimmi and Marmo, 2003, hereafter quoted as CM03). The larger ef- 
fects of asphericity are expected to occurr in homogeneous configurations, 
which have widely been investigated (Brosche et al., 1983; Caimmi et al., 
1984; CS90; Caimmi and Secco, 1992). Focaloidally striated density profiles 
involve far larger difficulty (e.g., Caimmi, 1995, 2003). 

The current investigation is mainly devoted to the following points: (i) 
expression of an equation of state for two-component systems; (ii) description 
of global properties deduced from selected density profiles; and (iii) applica- 
tion to elliptical galaxies belonging to a restricted sample, to be represented 
on the (O?/0) plane, for fiducial values of model parameters. 

The work is organized as follows. The basic theory of two-component 
systems with homeoidally striated density profiles, is reviewed and extended 
(to include both cored and cuspy matter distributions) in Section [21 The 
particularization to selected density profiles, involving explicit expressions, is 
made in Section[3l The results and related global properties are described and 
discussed in Section |H Elliptical galaxies belonging to a restricted sample. 
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are represented on the (Oj/0) plane, for fiducial values of model parameters, 
in Section where some considerations are drawn. The concluding remarks 
are reported in Section [6l 



2 Basic theory 

A general theory of two-component matter distributions has been exhaus- 
tively treated in an earlier paper (Caimmi and Secco, 1992), and the inter- 
ested reader is addressed therein and in parent investigations (MacMillan, 
1930, Chap. Ill, §76; Limber, 1959; Neutsch, 1979; Brosche et al, 1983; 
Caimmi et al., 1984) for further details. What is relevant for the current 
attempt, shall be reviewed and extended. If not otherwise stated, matter 
distributions should be conceived as continuous media instead of discrete 
particle sets (e.g.. Limber, 1959; Caimmi, 2007a). In the following, general 
definitions and related explicit expressions will be provided. Readers mainly 
interested to a simple reference case and a few cases of astrophysical interest, 
might directly go to Section [3], while readers mainly interested to the results 
and an application to elliptical galaxies, might directly go to Sections H] and 
[5], respectively. 

2.1 Kinetic-energy and self potential-energy tensors 

For an assigned density profile, the kinetic-energy tensor and the kinetic 
energy read (Binney and Tremaine, 1987, Chap. 4, §4.3): 



{Ei,in)pq = 7; Pi^i^^2,Xs)vpVgd^S ; p = l,2,3 ; g = l,2,3 ; (la) 



while the self potential-energy tensor and the self potential energy read 
(Chandrasekhar, 1969, Chap. 2, §10): 




s=l 



(lb) 




p= 1,2,3 



g = l,2 



,3 ; (2a) 




(2b) 



Js 
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where p is the density, d 5 = dxi (ix2 dxs an infinitesimal volume element 
placed at P = (xi,a;2,X3), VpVq = and vl = (v^) are arithmetic means 
calculated within d^S, Vpq and V are the gravitational tensor potential and 
potential, respectively (Chandrasekhar, 1969, Chap. 2, §10): 

Vpg(Xi,X2,X3) = G / p(Xi,X2,X3) ^ ^-3 dS ] 

•^^ R-W 
p= 1,2,3 ; g= 1,2,3 ; (3a) 

d^S' 



V{xi,X2,X3) = G p{x[,X2,X3)— — =^ ] (3b) 
■^s R-R' 

G being the constant of gravitation, R = OP and R' = OP', P' = {x'l, Xg, x'^), 
radius vectors with origin at the centre of inertia. 

Conformly to Eq. ([3]), the normalization of the potential satisfies the 
boundary condition to be null at infinite distances, and then positive else- 
where (e.g., MacMillan, 1930, Chap. II, §20; Chandrasekhar, 1969, Chap. 3, 
§17; Caimmi and Secco, 2003), instead of being null at the centre of inertia, 
and then negative elsewhere (e.g., Binney and Tremaine, 1987, Chap. 2, §1; 
Mouri and Taniguchi, 2003). 

Let i and j denote the subsystems of a two-component matter distribu- 
tion. The potential-energy tensor and the potential energy may be cast into 
the form (Caimmi and Secco, 1992): 

{EpQi)pg = [{Ei) sel]pq + [(yEij)int]pq + m.)int]pq+mUpq (4a) 
Epot = {Ei)scl + {Eij)int + {Eji)int + {Ej)sel ] (4b) 

where the expression of the potential energy (e.g., MacMillan, 1930, Chap. Ill, 
§76) has been extended to the potential-energy tensor; in addition, the in- 
teraction potential-energy tensor, [{Euv)mt]pq, and the interaction potential 
energy, {Euv)mt, read: 

[{Euv)int]pq = f Pu{Xl,X2,X3)[V^{Xi,X2,X3)]pqd^S ] 

u = i,j ; v = j,i ; (5a) 

(^™)int = -77 / pu{Xl,X2,X3)Vy{Xi,X2,X3)d^S ] 

2 JSu 

u = i,j ; v = j,i ; (5b) 

which are symmetric with respect to the exchange of one component with 
the other: 

[{Eij)mt]pq = [{Eji)int)]pq ; (6a) 

iEij)int = (-Eji)mt ; (6b) 
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for further details refer to earlier attempts (MacMillan, 1930, Chap. Ill, §76; 
Caimmi and Secco, 1992). 

The tidal potential-energy tensor and the tidal potential energy read (e.g., 
Caimmi and Secco, 1992): 



f dVv 

[{Euv)tid]pq = Pu{Xi,X2,X3)Xp-^d'^S ; 



u^ij ; v^j,i ; (7a) 
{Euv)tid= Puixi,X2,X3)J2^s-K-^<i!^S ; 

^ Su S=l 

u^ij ; v^j,i ; (7b) 

where the tidal potential energy, (£'uv)tid, may be conceived as the virial of 
the M-th component in connection with the tidal field induced by the I'-th 
component (Broschc ct al., 1983). 

The tensor and the scalar virial theorem for a single subsystem, within 
the tidal field induced by the other one, read (Caimmi et al., 1984; Caimmi 
and Secco, 1992): 

2[(-E„)kin]pQ + [(-E„)sel]pg + [(-E„t,)tid]pg = ; u = i,j ] v=j,i ] (8a) 

2(£^„)kin + (^«)sei + (^«,;)tid = ; u ^ ij ] v^j,i ] (8b) 

which is the generalization of previous results related to one-component sys- 
tems (e.g., Chandrasekhar, 1969, Chap. II, §11; Binney and Tremaine, 1987, 
Chap. 4, §3). The vahdity of the relations (Caimmi and Secco, 1992): 

[{Eij)tid]pq + [(-£'ii)tid]pg = [{Eij)i^t]pq + [{Eji)ij^t)]pq ', (9a) 
{Eij)tid + {Eji)tid = {Eij)int + {Eji)mt ; (9b) 

implies the following: 

[(-^ut;)tid]pg [(-^ui))int]p5 ~l~ [(-^uti)res]pg j 

u^i,j; V ^ j,i ; (10a) 
{Euv)m ^ {Euv)int + {Euv)res u^ij ; v^j,i ; (10b) 

where [{Euv)res]pq and (i?„„)res are the residual potential-energy tensor and the 
residual potential energy, respectively, which are antisymmetric with regard 
to the exchange of one component with the other: 

[(-S'ij )rcs]pg [(-^j'i )rcs]p(; > (11a) 

(-Ejj)res = —(Eji)rcs (H^) 

for further details refer to earlier attempts (Caimmi and Secco, 1992; Caimmi, 
2007b). 
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It can be seen from Eqs. ([H]), ffTUl) . and ffTT]) . that the tidal potential- 
energy tensor, [{Euv)tid]pq, and the tidal potential energy, {Euv)tid, are made 
of a symmetric term, [{Euv)int]pq and {Euv)mt, and an antisymmetric term, 
[{Euv)res]pq and (-E'„„)res, respectively, with regard to the exchange of one 
component with the other. 

The virial theorem in tensor and in scalar form, expressed by Eqs. ((H]), 
may be cast into the more compact form: 

2[(^n)kin]pg + [(^™)vir]p<? = ; u = i, j ; y = j,i ; (12a) 
2(-E„)kin + (-E™)vir = ; u = i,j ] v=j,i ] (12b) 

where the virial potential-energy tensor, [{Euv)vir]pq, and the virial potential 
energy, {Euv)vir, are defined as: 

[{Euv)viT]pq = [{Eu)scl]pq + [{Euv)tid]pq 'i U = i, j ; V = j,i ] (13a) 

{Euv)vh = {Eu)sci + {Euv)tid ; u = i,j ; v = j,i ; (13b) 

where, in general, the virial potential energy is usually named "the virial of 
the system" (Clausius, 1870). In the case under discussion, the "system" 
relates to the u-th component within the tidal potential induced by the v-th 
component. 

For assigned density profiles, the virial potential-energy tensor and the 
virial energy of each subsystem can be determined, together with their kinetic 
counterparts via Eqs. fll2al) and (112b[l . respectively, which constrain, in turn, 
the orbital anisotropy on each subsystem. More specifically, the sum of mean 
orbital kinetic-energy tensor and kinetic energy, has to reproduce [(-En)kin]p(j 
and (-Eu)kin, respectively. 

When a system is not entirely included within its truncation radius, the 
usual form of the virial theorem, 2_Ekin + Epot = 0, should be extended 
as 2£'kin + i^pot = 3pS, where the last is a surface term (e.g.. The and 
White, 1986; Carlberg et al., 1996; Girardi et al., 1998). In dealing with 
a specified subsystem, it shall be intended that no mass exists outside the 
related truncation radius, which makes a null surface term. 

To avoid the determination of the gravitational potential, which is the 
most difficult step towards an explicit expression of the potential-energy ten- 
sors and potential energies, a particular procedure shall be followed under 
the restrictive assumption of homeoidally striated density profiles (Roberts, 
1962). 
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2.2 Homeoidally striated density profiles 

Let the isopycnic (i.e. constant density) surfaces be defined by the following 
law (CMOS): 



p = pV(0 ; /(i) = i ; (14a) 

e = j:^, ; 0<e<S ; (14b) 

where = p(l) the density and the semiaxes, respectively, of a 

reference isopycnic surface, and S corresponds to the truncation isopycnic 
surface, related to semiaxes, a^. The scaled radial coordinate, ^, and the 
scaled density, /, may be conceived as the generalization of their counterparts 
related to polytropes (e.g., Chandrasekhar, 1939, Chap. IV, §4; Horedt, 2004, 
Chap. 2, §2.1). 

According to Eqs. (|T^ . the scaling density, p^, and the scaling radius, r\ 
correspond to a single boundary, which allows the description of both cored 
and cuspy density profiles. The assumption that the system is homeoidally 
striated implies the relation (CMOS): 



and, in particular: 

XL . . 

which, in any case, is independent of the radial coordinate, r, of the generic 
point on the selected isopycnic surface. 

2.3 Mass and inertia tensor 

The function (Roberts, 1962): 

F{0 = 2 r fiO^'d^' ; (17a) 



5 

F{^) = ; ^ = -2^/(0 ; (17b) 

r fiOC = ^ r FiOC-' ; n>l ■ (17c) 
JO 2 Jo 

allows the calculation of the total mass as (CM03): 

M = z/^,sM^ ; (18a) 
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^ ^V(Ode ; (18b) 



2 Jo 

Mt = ^pV,44 ; (18c) 

and the inertia tensor as (CMOS): 

Ipg = 5pgiyi^rM^{alf ; (19a) 

^inr=^ pFiOed^ ; (19b) 

2 JO 

where the coefficients, i^mas and z/jnr, are shape-independent and may be con- 
ceived as proffie factors (Caimmi, 1993), and 6pq is the Kronecker symbol. 
For a different formulation of the inertia tensor refer to earlier attempts (e.g.. 
Landau and Lifchitz, 1966, Chap. VI, §32; Bett et al., 2007). 

The mass within an isopycnic surface, p = p"^ f{^), is M(^) = z/mas(0^^ 
where z/mas(0 is expressed by Eq. (118bll with ^ instead of S. The related 
mean density is p{C,) = M{C,)/S{C,), where S{^) is the volume bounded by 
the isopycnic surface. 

2.4 Potential-energy tensors and potential energies 

The self potential-energy tensor and the self potential energy read (CM03): 

G(Mt)2 

{Esel)pq = —SpgUsel 1 Bp ] (20a) 

G(Mt)2 

E,,i = . B ■ (20b) 

i^.ei = ^ r^'(Ode ; (20c) 



16 Jo 



Bp = ep,eps / (1 + .')-^/^(l + + els')-'^' ds' ; (20d) 

3 a 
B = Y,Bp ■ e^n = — \ (20e) 



s=l 



where ^sei is a profile factor, emn are axis ratios, and Bp are profile factors 
which, ipso facto, depend on the axis ratios only (Caimmi, 1992). 

The above results are related to a single subsystem: strictly speaking, all 
the quantities defined in the current section should be labelled by the index, 
M, in connection with the M-th subsystem (e.g., Caimmi and Secco, 1992), 
but it has been omitted to gain clarity. On the other hand, the formulation 
of other potential-energy tensors and potential energies necessarily involves 
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(at least) two components. The related calculations are very difficult in the 
general case, and for this reason only the special situation of similar and 
similarly placed boundaries shall be considered. 

With this restriction, the combination of Eqs. ( fT5i) . ( flGil . and f llScp . re- 
lated to both the density proffies, the inner to be denoted as i and the outer 
as j, yields: 

^^ = y% ; =^ = ^ ; Mii^ = ^ ; (2ia) 

y = : 'u' = Hr: m = —7- : m' = — ^ : (21b) 

which allows the expression of the other potential-energy tensors and poten- 
tial energies, as (CMOS): 

G(Mt)2 

[(-£'tiv)xxx]p(; (^pq ~ +~ {yuv)xKxBp , (22a) 

(-^™)xxx z +~ ('^Mt))xxx-S 1 (22b) 

u = i,j ; V = j,i ; xxx = int, tid, res, vir ; (22c) 
and the explicit expression of the profile factors reads (CM03): 
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w;('°*)(^) ; (23a 



(^..)tid = -^mt«;(-*)(r^) ; = -l^w^'^'^Hv) ; (23b) 

O o 771' 

{i^ijU = -i^jiU = -^mt [w^'-'\ri) - w^'^^'\ri)] ; (23c) 
(i^™)vir = («^n)soi + («^Mt,)tid ; u = i,j ] v = j,i ] (23d) 
V = ^ = ^ ; (23e) 

y 

where the functions, w^"^*) and w^^^^\ are defined as (CM03): 

/V,(0)^OdO ; (24a) 
rv d F 

w^^^'\v)= / m^)-7f^,d^J ; (24b) 

Jo (iti 



in conclusion, Eqs. fl20p - (P^ allow the calculation of the potential-energy 
tensors and potential energies for homeoidally striated density profiles related 
to similar and similarly placed boundaries. 
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In the case under discussion, due to Eqs. f l22|) . the fractional virial potential- 
energy tensor component equals the fractional virial potential energy, as: 

_ [(-^j»)vir]pg _ {Ejijvh _ {fn^) (^ji)vir _ (25) 

which, for assigned density profiles, depends on the reference fractional mass, 
, and the fractional scaling radius, y^, according to Eqs. fl2ip - (P^ . The 
substitution of Eqs. (!?!]) into fl^ yields: 



y Si 



[jyi 



i^jihr . ^26) 



] J mas J V * J y vir 



which, for assigned density profiles, depends on the fractional mass, m, and 
the fractional truncation radius, y, according to Eqs. ( 12T]) -(12^. 

Strictly speaking, Eqs. ( l2Tl) -( l26l) are valid provided the indices, i and j, 
denote the embedded and the embedding subsystem, respectively, which im- 
plies y > I. This is why the above procedure is valid only for the inner com- 
ponent, where the Poisson equation instead of the Laplace equation holds 
with respect to the gravitational potential. The remaining results, related to 
the outer component, are obtained by use of the symmetry of the interaction 
potential-energy tensor and interaction potential energy, Eqs. ([6]), and the 
antisymmetry of the residual potential-energy tensor and residual potential 
energy, Eqs. ffTT]) . 

If the role of the two subsystems is reversed, < ?/ < 1, it has to be kept in 
mind that, in this case, the inner and the outer component are denoted by the 
indices, j and i, respectively, and the reversion must be done in Eqs. ( l2Ti) -( l26l) . 
through the following steps: (i) make the changes: m m~^; m) (m"'')"^; 

y y'^; y^ {y^y^; ^ Ej-, (ui)^^ ^ (i^jOmas; ('^i)sei ^ ('^jOsei; (h) 

for assigned y > 1 and y"^ = {Ei/Ej)y, calculate the functions, w^^^^\ri) 
and w^'^'^^^T]); (iii) calculate the profile factors, (t'jj)tid and {h'ji)tid, and the 
remaining ones when needed; (iv) calculate the fractional virial potential 
energy, (j); (v) make the changes, m — > m~^; m'^ (m"'')^^; y ^ y^^; y'^ 
(?/^)~^; (f) —>■ (f)^^] which allow the extension of the fractional virial potential 
energy, (f) = {Eji)^^^/ {Eij)^-„, to the domain, < ?/ < 1. 

In absence of truncation radius, S — >• -|-oo, t] +oo, the reversion occurs 
when the density drops to zero and nothing changes except in infinitesimal 
terms of higher order and infinite terms of lower order. Accordingly, there is 
no need to perform the reversion in this case. 

The relation, = (p{y\ m^), or its counterpart, = (f){y, m), expressed by 
Eqs. ( |25ll and (l26l) . respectively, may be conceived as an equation of state for 
two-component systems with assigned homeoidally striated density profiles. 
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At the price of a major complexity, it can be formulated for any kind of 
two-component systems, hereafter quoted as "two- component macrogases", 
or more shortly as "macrogases" . 

3 Special cases 

The explicit expression of the macrogases equation of state is, in general, 
rather cumbersome and numerical computations should be preferred to this 
aim. On the other hand, the related procedure is conceptually simple, as the 
functions, w^^^^\r]) and w^^^^\r]), and the fractional virial potential energy, 
4>, may be calculated using Eqs. (!2U1) - (PUI) . 

Aiming to provide a description of the general trend and related features, 
a limited number of simple possibilities shall be analysed in detail, selecting 
density profiles from the family: 

•^^^) = ^7(l + ^")X ' ^ = ' ^^^^ 

which is defined by three parameters, 7). For further details refer to 

earlier attempts (e.g., Hernquist, 1990; Zhao, 1996; Caimmi and Marmo, 
2004; Caimmi et al., 2005; Caimmi, 2006b), where special cases were fully 
investigated. 

In the following, the macrogases equation of state shall be determined for 
a simple but unrealistic density profile, to be taken as a reference case, and 
a few density profiles of astrophysical interest. The reader whose attention 
is mainly directed to the results and/or the astrophysical applications of the 
model, is free to jump directly to Section H] and/or El respectively. 

3.1 UU macrogases 

The related density profiles maintain uniform, (a, /5, 7) = (0, 0, 0), and Eq. fl27|) 
reduces to: 

A(a) = i ; 0<e«<H„ ; u = i,j ; (28) 

which is equivalent to polytropes with index, n = (e.g., Chandrasekhar, 
1939, Chap. IV, §4; Caimmi, 1986) but implies unphysical situations for stel- 
lar fiuids (Vandervoort, 1980). The particularization of Eqs. f ll7ap . ( ]18b[) . 
f l20cl) . f l24al) . and ( ]24bl) to the case of interest yields: 

Fu{U = El-eu ; u = z,j ■ (29) 
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— ^3 . 

ujmas • — 'u 1 

) sel 



10"" ' 
4 



u = i,j ; 
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(30) 
(31) 

(32) 

(33) 



and, using Eqs. ([230) and f l23ep . the UU macrogases equation of state follows 
from the particularization of Eq. (!26|) to the case of interest, as: 



(m 



t^2 



1 + i^lfV- 

(mt)y5 ^2^ 



1 + 



2/> 1 



(34) 



and the extension of the above function to the domain, < y < 1, following 
the procedure outlined in Subsection 12.41 yields: 



1 + 



5 3 2 
2 ~ 2^ 



0<y <1 



(35) 



where it can be seen that in the special case, y = 1, Eqs. ( IMl) and ([35 
do coincide. If, in addition, Sj = Ej, which implies 
Eqs. ([IHD and ([HD, then Eqs. and ([SSD reduce to: 



y, m 



m, via 



m = m 



y = y 



1 ; 



(36) 



or 0(1, m) = m. 

Owing to Eqs. (1341) and (I35p . the fractional virial potential energy, 0, is 
independent of the fractional truncation radii, (S,, Ej). 



3.2 PP macrogases 

The related density profiles (Schuster, 1883; Plummer, 1911) imply (a, /5, 7) = 
(2, 5, 0), and Eq. ([27]) reduces to: 

25/2 

fui^u) = _^ ^2)5/2 ' 0<^„<S„ ; u = i,j ] (37) 
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which is equivalent to polytropes with index, n = 5; for a formal demonstra- 
tion, see Appendix |X1 For one-component systems, hydrostatic equilibrium 
holds regardless from the nature of the fluid (Vandervoort, 1980). 

The particularization of Eqs. (ITTall . (flSbll . (l20cll . (l2lall . and fl24bl) to the 
case of interest yields: 

27/2 



( ) mas 



25/2'^3 



, , S^3S^^8H|+13) „ 
l^^ujsei = , ^2^3 h 3arctan^„ 



u = i,j ; (38) 
(39) 

u = i,j ; (40) 



^(mt)(^) = - — (yt)2 1 [(^^^ + 7my\ ia) + 4[(yt)2 _ ly.pjy^ ia) 



+ 



3[(l/t)2_l]3 
1 



3[(l/t)2 _ l]3(^2 + l)[(yt)2,^2 + 1]2 (1 + ^2)3/2 3[(?/t)2,^2 + l]3/2 

pM)(r^,yt) = r7w'r^2 + i /(^t)2^2 + i 



X [{y^fv'iv' + 1) + (y^)'r]'(rr/^ - 4) + (i/t)2(ii^2 ^ 5) ^ 3] . 
a = arcsinh?] ; 7^ 1 ; 

128 



(mt)(^) 



t^2 



arctanr/ ^ r]{37]'^ + 877^ — 3) 



16 



1 



(l + S2)3/2 3(,^2 + 1)3/2 



2/ 



48(7^2 + 1)3 

^-1 ; 



;(41a) 



(41b) 
(41c) 



(41d) 



128 r [7(i/t)2 + i]ii?(^t, i«) + [3(yt)4 _ 2(^t)2 _ i]ip(^t^ i^) 



3[(l/t)2-l]3 



3[(yt)2 _ 1]3(^2 + l)2[(^t)2^2 + 1] (1 + ^2)3/2 3(^2 + 1)3/2 

X [(yt)4(7r^4 ^ 11^2 ^ 3) ^ (^t)2(^4 ^ 4^2 ^ 5) ^ ^2] . 
a = arcsinh?7 ; 7^ 1 ; 

^(ext)(^)^_l^ 



arctani] ^ ■r]{3ri + 8r/ — 3) 



16 



(1 + S2)3/2 3(,^2 + 1)3/2 



48(7^2 + 1)3 

2/^ = 1 ; 



(42a) 



(42b) 
(42c) 



(42d) 
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where i is the imaginary unit, F and E are incomplete elhptic integrals of 
the first and second kind, respectively, defined as: 

Fik,p)^r .. - r '1 ; (43a) 



Vl - sin" e Jo VT^Vl - kH'^ 

rf3 I rx ,/l _ 1.2+2 

E{k,(3)= Jl-k^sm^ede= ^ , dt ; (43b) 

JO Jo vl — t"^ 

X = sin /3 ; t = sin 6* ; (43c) 

for further details refer to specialized textbooks (e.g., Spiegel, 1968, Chap. 4, 
§§34.1-4). 

Using Eqs. (123bp and the PP macrogases equation of state is 

obtained from the particularization of Eq. ( !26l) to the case of interest for the 
domain, y > 1. The extension to the domain, < y < 1, can be done 
following the procedure outlined in Subsection 12. 4[ 

In absence of truncation radius, the density drops to zero when the radius 
goes to infinity, H +oo, rj +oo, and Eqs. flHSl) - fH21) reduce to: 



27/2 I 

=«^+oo 3(1 + iiyi^ 

_ lim {uu)mas = 2^/" ; u = i,j ] (45) 
Sn 

^lim (i/„)ggi = — ; u = i,j ; (46) 



lim w^'"^*) 

r;— >+oo 



iv) 



128 t.2 [jy^f + 7]E{k, ir/2) - [5{y^f + 3]F(fc, 7r/2) 
3 ^ 3[(yt)2 _ i]3 

(y¥ ; 1 ; (47a) 



lim w('"*)(77) = 



t 



(47b) 



lim . 128 t + 5]F(fc, vr/2) - [7{y^)' + l]E{k, 7r/2) 

r;-+oo ^"l^ 3 ^ 3[(yt)2_l]3 

A; = yr^I7(^ ; _^ 1 . (48a) 

Ajr 

lim w('="*'(77) = ; y^ = 1 ; (48b) 

r;^+oo 3 

where, in particular, the related expression of the tidal energy coincides with 
its counterpart calculated in an earlier attempt (Valentinuzzi, 2006, Chap. 4, 
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§4.1). Using Eqs. fl23bl) and fH3|) - fHH|) . the PP macrogases equation of state in 
the special situation under discussion, is obtained from the particularization 
of Eq. ( !26l) to the case of interest, for the domain, y > 0. 

If, in addition, y = y'^ = 1, Ej = Ei, the combination of Eqs. (l26il and 
IHJ-dlHD yields Eq. (I36l). 



3.3 HH macrogases 

The related density profiles (Hernquist, 1990) imply 7) = (1,4, 1), and 
Eq. ( 1271) reduces to: 



< C„ < S„ ; u = i,j 



(49) 



which has been proved to be consistent with nonnegative distribution func- 
tions, in the parameter range of interest (Ciotti, 1996). 

The particularization of Eqs. (fTTal) . ffTSbl) . (IM . (Ma^ . and fl24bD to the 
case of interest yields: 
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;i + e„)2 (1 + HJ2 ' 

u = i,j ; 



u = t,J 



12-2 



(1 + 



(50) 
(51) 
(52) 



^2(t/t-l)4 
yirj + 1 yTr^ -f 1 



(yt - l)2^t^ ^ 2(yt - l)r^ 



1 + r/ 



1 



2{1 + E,y{y^y 



1 - 



U7('°*)(^) = 

1 1 



-128 



12 



2y^r] + 1 
4r7 + 1 



+ 1 



T] 



2(1 + S,)2(r^ + 1)^ 



(7^ + 1)4 

; 1/^ = 1 ; 



(53a) 



(53b) 



-128 



2{y^~iy 



(l/t - l)2r^ ^ 2(yt)2(^t _ 1)^ 



(^ + 1) 



1 + y'^r] 



16 



ri + l 



2(1 + S,)2(r7 + 1)^ 



w 



(cxt) 



128 <^ — 
I 12 



4r/ + 1 
(^ + 1)' 



- 1 



(54a) 



using Eqs. ( 123bl) and (!5Tl) -( l5^ . the HH macrogases equation of state is ob- 
tained from the particularization of Eq. fl2U|) to the case of interest for the 
domain, y > 1. The extension to the domain, < y < 1, can be done 
following the procedure outlined in Subsection I2.4[ 

In absence of truncation radius, the density drops to zero when the radius 
goes to infinity, H — > +oo, r] — > +oo, and Eqs. (I50 i) -( 1M1) reduce to: 



_ lim Fui^u) 



lim fz/„ 



(l+^n)^ ' 

12 ; u 



u 



lim (u, 



u Jsel 



12 



u 



= hJ 



(55) 
(56) 
(57) 



hm w(^'^*)(r7) 
lim w^'^^^Hn) 



64yt 



(yt - 1)^ 
32 



-2(2|/t + i)lnyt + (^t_i)(^t + 5) 



(58a) 
(58b) 



lim w(^"*)(r7) 
lim w("^*H^) 



64 



(y^ 
32 



2yt(^t + 2) lnyt_(^t_i)(5^t + i) 



(59a) 
(59b) 



where, in particular, the related expression of the tidal energy coincides with 
its counterpart calculated in an earlier attempt (Valentinuzzi, 2006, Chap. 4, 
§4.2.1). Using Eqs. (123bp and (!56|) - (!59|) . the HH macrogases equation of state 
in the special situation under discussion, is obtained from the particulariza- 
tion of Eq. 0261) to the case of interest, for the domain, y > 0. If, in partic- 
ular, y = = 1^ 5^. 

Eq.(l36l). 



Sj, the combination of Eqs. ( ]26l) and (I56l) -(l59l) yields 
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3.4 HP macrogases 



The inner density profile (Hernquist, 1990) implies 7) = (1,4, 1), which 
is defined by Eq. ( l49i) . and related functions and parameters by Eqs. (!50l)-(l52l) 
or ( l55l) -( l57l) in the special case of no truncation radius. 

The outer density profile (Schuster, 1883; Plummer, 1911) implies (a, [3, 7) 
(2,5,0), which is defined by Eq. (1H7I) . and related functions and parameters 
by Eqs. (!HS|) - P(1 or (H^ - fHU|) in the special case of no truncation radius. 

With regard to the tidal potential-energy terms, the particularization of 
Eqs. fl24ap and (124bl) to the case of interest yields: 



w 



(int) 



$r)(r^) 



128^2 



{y 



t^4 



(int) 



2[(yt)2 + l]3 

3yt[3(yt)2_2] 



2 

- + 



(int) 



- 2 (in 

p(int)j^^-j 



,)3/2 



In 



2[(yt)2 + l]3(^t,^+ 1)2(^2 + 1)1/2 



2[(yt)2 + l]7/2 ^^^^ ^ ^ y'(yt)2 + ij ' 

-2(t/t)5^(2772 + 1) + (^t)4(^2 _ ^) ^ (yt)3^(11^2 ^ 

-4(yt)2(4^2 + 3) + 2|/t^-2 ; 



(60a) 
(60b) 

(60c) 
(60d) 



w 



(cxt) 



-64v^ 
ll(l/t)3 



1 



1 



+ 



3(1 + H, 

p(ext) j^^) 



1+772)3/2 



3[(l/t)2 + l]3 3[(l/t)2 + l]3(^tr7 + l)(^2 + 1)3/2 



(yt)2[2(^t)2_3] 



In 



(61a) 
(61b) 

(61c) 



[(yt)2 + l]7/2 ^^^^ ^ ^ ^(yt)2 + 1 

(yt)y (2r/2 + 3) - (yt)^r7(4r/2 + 5) - [y^fiUr]" + 217]^ + n) 
-{y^fr]{'iTf + h)+y\r]^ + QTf + A)+Tf ; (61d) 



using Eqs. fl23b1) . (1391). fl40|) . ([5lD, and ([52]), the HP macrogases equation of 
state is obtained from the particularization of Eq. (1261) to the case of interest 
for the domain, y >1. 

In absence of truncation radius, the density drops to zero when the radius 
goes to infinity, H -> +cx), t] +oo, and Eqs. ([38]) -(SO]), ([50]) -([521), reduce to 
([4l])-([46]), ([55])-([57D, respectively, and Eqs. ([nn])-([6l]) reduce to: 



lim w;('°*)(r7) 



QA^y 



3[(yt)2 + ip 



- A{y^f - I2{y^f + lly^ + 2] 



,t^3 



t^2 



,t 
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3yt[3(yt)2_2] 



In- 



1 - 1 



(yt)2 + 1 



lim w('°*)(r7) = - 



+ ^(y^y + 1 
72-1' 



2^2 + 3 In 



1/^ = 1 ; 



(62a) 
(62b) 



lim w(^"*)(r7) 



64^2 



3[(yt)2+l]3 



[2(yt)4 + ii(^t)3_i2(^t)2_4^t + i] 



t^2 



3(yt)2[ -2(yt)2 + 3] + 1-1 



lim w('"*)(r7) 



(t/t)2 + 1 

c r 

2V2 + 31n 



y^ + ^{y^r + 1 

v^- 1 



V2 + 1 



(63a) 



(63b) 



where, in particular, the related expression of the tidal energy coincides with 
its counterpart calculated in an earlier attempt (Valentinuzzi, 2006, Chap. 4, 
§4.2.2). Using Eqs. fl23bl) and , 1^ , 1^ , 1^ , §^ , ^ , the RF macro- 
gases equation of state in the special situation under discussion, is obtained 
from the particularization of Eq. ( |26l) to the case of interest, for the domain, 

y>o. 



3.5 HN macrogases 

The inner density profile (Hernquist, 1990) implies (a,/?, 7) = (1,4, 1), which 
is defined by Eq. (149!) . and related functions and parameters by Eqs. fl50|) - (!52|) 
or (!^ - (!Fn) in the special case of no truncation radius. 

The outer density profile (Navarro et al., 1995, 1996, 1997) implies (a, /?, 7) = 
(1, 3, 1), and Eq. (1271) reduces to: 

= |;(r|yi : i (64) 

which, together with its H counterpart expressed by Eq. (HUI) . have been 
proved to be consistent with nonnegative distribution functions, in the pa- 
rameter range of interest (Lowenstein and White, 1999). 

The particularization of Eqs. f ll7al) . ( 118bl) . and fl20cl) to the case of interest 
yields: 
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ln(l + Ej 



1 + H 



'J J 



sel 



36 



S,(2 + H,)-2(1 + H,)ln(l + H, 
(1 + 2,)^ 



(66) 
(67) 



with regard to the tidal potential-energy terms, the particularization of Eqs. (I24al) 
and fl24bp to the case of interest yields: 



+2 In 



6%t 



(2/t-l)2^t^(^t^ + 2) , 2ytry(yt_l) 



(yt - 1)3 
r] + l (y^-l) 



yt (i/tr; + 1)2 



?/1"r7 + 1 



(1 + S.)^ 
wM)(ry) = 



64r/2 



V + 3 1 

3(77+1) 1 + 



64 



y'^v 



; y^ = i 



(68a) 
(68b) 



{y'f — ly [ r] + 1 y'^rj + 1 — 1 y'^rj + 1 



+ ln(r/ + 1) 



-64 



7] + 1 

'lri^{ri + 3) 1 
6 (77 + 1)3 ~ (1 + S,)2 



7/+ 1 



+ ln(7/ + 1) 



y 



t-1 



(69a) 



(69b) 



using Eqs. fl23bp . ([51]), fl3^ . flUBl) . and flUTj) . the HN macrogases equation of 
state is obtained from the particularization of Eq. 0261) to the case of interest 
for the domain, y > 1. 

In absence of truncation radius, the density drops to zero when the radius 
goes to infinity, H — > +oo, t] +oo, and Eqs. (I50l)-(l52l). reduce to ( l55l) -( l57l) . 
respectively, and Eqs. reduce to: 



(70) 
(71) 
(72) 



lim (z/j)mas = +00 ; 



lim (z/, 



S,-— >+oo 



j J sel 



36 



lim m;('"^*)(77) 
lim m;("^*)(^) 



64 



(7/t - 1)3 

64 



1 ; 



7/t^l ;(73a) 
(73b) 
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lim w^'''''\7]) 



64 



lim = - — 

3 



?7— *+oo 



_ 1 



; (74a) 
(74b) 



where the self potential-energy profile factor remains finite, although the 
mass profile factor undergoes a logarithmic divergence. Using Eqs. (I23bl) and 
(!72l) - (!69l) . the HN macrogases equation of state (where mJ and appear 
instead of Mj and m) in the situation under consideration, is obtained from 
the particularization of Eq. 0261) to the case of interest, for the domain, y > 0. 



4 Results 

The macrogases equation of state, Eq. (1261) . is represented in the next figures 
for a number of cases analysed in Section [31 In particular, owing to Eqs. flMj) 
and fl3^ . the UU macrogases equation of state is independent of the scaled 
truncation radii, (Hj,Hj). 

The macrogases equation of state in absence of truncation radius, is plot- 
ted in Fig.[I] for cases (from top left in clockwise sense) UU, PP, HH, HP, 
where = Sj and (from bottom to top in each panel) m =1, 2, 6, for 
cases UU, PP, HP, and m =10, 20, 60, for case HH. Shallower density 
profiles (UU, PP) show the occurrence of two extremum points: one mini- 
mum on the left and one maximum on the right, for arbitrary values of the 
fractional mass, m. Steeper density profiles (HP, HH) still exhibit extremum 
points, but none below a treshold, where the critical isofractional mass curve 
shows a single horizontal inflexion point. 

As already noticed in an earlier attempt (CS90), the same trend is shown 
by van der Waals' (1873) isothermal curves and, in fact, the van der Waals' 
equation of state looks similar to the macrogases equation of state, Eq. ( l26i) . 
where three variables also appear. Thick curves represent the locus of mini- 
mum (left branch) and maximum (right branch) points. The critical isofrac- 
tional mass curve, when it appears, is also thickened. 

A main feature is that, above the critical isofractional mass curve, for any 
selected value of m, a range in fractional virial potential energy, 0, exists, 
0min < (p < 0max5 where three different configurations i.e. different fractional 
truncation radius, y, correspond to the same value of 0. By analogy with van 
der Waals' isothermal curves, it could be argued the existence of a bell-shaped 
region in the Oycp plane, within which the extremum points are located, and 
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where a phase transition occurs. Further investigation is needed on this 
point, but it hes outside the aim of the current paper. 

The macrogases equation of state in presence of truncation radius, is 
plotted in Figs. [2] and [31 respectively, for different choices of the fractional 
truncation radii, represented on each panel where (from bottom to 

top) m =10, 20, 60, and approximate values of the parameters related to 
the critical point i.e. the horizontal inflexion point on the critical isofractional 
mass curve (marked by a St. Andrew's cross), are listed in Tab.[Tl 

It can be seen that the occurrence of the truncation radius makes little 
change to the trend of the isofractional mass curves. In general, decreasing 
the outer scaled truncation radius, Ej, yields more pronounced extremum 
points and vice versa. Much smaller changes, in the same sense, are produced 
by decreasing the inner scaled truncation radius, Hj. This is why, in the cases 
under discussion, the mass of the outer subsystem is dominant with respect to 
the mass of the inner subsystem, m 3> 1. The critical isofractional mass curve 
does not appear for earlier truncated HN density profiles, (Sj,Sj) = (5,5), 
which are shallower in comparison with their latter truncated counterparts. 

Reduced truncation radii, S = R/r\ may be conceived as concentrations 
of the related matter distribution (e.g., Navarro et al., 1997; CMOS; Caimmi 
et al., 2005) and the selected range, 5 < H < 20, is consistent with the 
results of dark matter (hereafter quoted as DM) halo numerical simulations 
(e.g.. Bullock et al., 2001) and elliptical galaxy (hereafter quoted as EG) 
observations (e.g., Lowenstein and White, 1999). 

5 Application to elliptical galaxies and their 
hosting haloes 

According to current cosmological scenarios (e.g., Mota and van de Bruck, 
2004; Percival, 2005; Horellou and Berge, 2005; Maor and Lahav, 2005; Nunes 
and Mota, 2006), density perturbations at recombination epoch {z ^ 1100) 
initially expand with the universe, turn around, collapse, and finally virialize 
(at least in their inner and denser regions). 

Virialized density perturbations, such as typical elliptical galaxies (EGs) 
and clusters of galaxies, may be idealized, to a first extent, as two homeoidally 
striated, similar and similarly placed, density profiles. In the following, at- 
tention shall be focused on EGs. 
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Table 1: Approximate values of parameters related to the critical point i.e. 
the horizontal inflexion point on the critical isofractional mass curve, for the 
density profiles under investigation. Bars instead of numbers mean absence 
of critical curve in the corresponding case. Numbers in brackets denote the 
values of scaled truncation radii, (Ei.Ej). In absence of truncation radius, 
S — +00, the case considered is r.j/Si 1, which implies = y. For 
homogeneous configurations (case UU) the isofractional mass curves are in- 
dependent of 
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5.1 General considerations and main assumptions 

A recent investigation performed on an optically complete sample of 42 EGs, 
for which X-ray gas temperatures and central stellar velocity dispersions 
were determined (Davis & White 1996), has shown evidence that, in general, 
EGs contain substantial amounts of DM (Loewcnstein and White, 1999). 
Accordingly, more than about 20% and 39%-85% of the total mass within 
one and six optical radii, respectively, is in form of (non baryonic) DM, 
depending on the stellar density profile and observed value of X-ray gas 
temperature and central stellar velocity dispersion. The comparison between 
the velocity dispersion distributions for DM and stars, assuming isotropic 
orbits, shows that the DM is dynamically "hotter" than the stars, by a factor 
1.4-2 (Loewenstein and White, 1999). 

The above investigation cannot be considered as conclusive in favour of 
the existence of DM haloes hosting EGs. In fact, it has been pointed out that 
the attenuation (in particular, the scattering) by dust grains has the same 
effect on the stellar kinematics as a DM halo (Baes and Dejonghe, 2001). 
According to a recent attempt, no strong evidence for DM haloes within 1-3 
optical radii has been found in a restricted sample of 4 early-type EGs, using 
dynamical modelling (Samurovic and Danziger, 1995). Beyond 1-3 optical 
radii, the X-ray methodology shows the need of DM where an X-ray halo is 
detected (Samurovic and Danziger, 1995). A similar trend is exhibited by 
an additional early- type EG (Samurovic and Danziger, 1996). In any case, 
current cosmological scenarios (ACDM, QCDM) predict DM haloes hosting 
EGs, as well as spiral galaxies, for which there are lots of empirical evidence 
(e.g., flat rotation curves well outside optical radii). For this reason, EGs are 
also assumed to be embedded within DM haloes. 

An analysis on the evolution of the physical properties of cosmological 
baryons at low redshifts {z ^ 5) has recently been performed (Valageas et 
al., 2002), which (i) yields robust model-independent results that agree with 
numerical simulations; (ii) recovers the fraction of matter within different 
phases and the spatial clustering computed by numerical simulations; (iii) 
predicts a soft X-ray background due to the "warm" intcrgalactic medium 
component, that is consistent with observations. The related baryon fraction 
in the present universe is found to be 7% in hot gas, 24% in the warm inter- 
galactic medium, 38% in the cool intergalactic medium, 9% within star-like 
objects and, as a still unobserved component, 22% of dark baryons associated 
with collapsed structures, with a relative uncertainty no larger than 30% on 
these numbers. Then the amount of still undetected baryons is about one 
flfth of the total, one fourth of the observed baryons (intergalactic medium, 
stellar components, and hot gas), and at least twice the stellar-like compo- 
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nent. 

According to recent investigations, HH (e.g., Ciotti, 1996) or HN (e.g., 
Loewenstein and White, 1999) models provide viable representations for the 
inner, baryonic, and the outer, nonbaryonic subsystem, respectively. The 
above mentioned mass distributions were found to be self-consistent, in the 
parameter range of interest, with regard to the non negativity of the distri- 
bution function (e.g., Ciotti, 1996; Loewenstein and White, 1999) by use of 
a theorem stated in an earlier attempt (Ciotti and Pellegrini, 1992). 

If undetected baryons in EGs are present as hot gas, the gaseous sub- 
system is expected to be less concentrated than the stellar one, as in the 
Coma cluster of galaxies (e.g., Briel et al. 1992). If undetected baryons in 
EGs are present as unseen stars, the undetected subsystem is expected to be 
as concentrated as the stellar one. Either assumption is necessary in deal- 
ing with two-component systems. The real situation lies between the above 
mentioned limiting cases. 

The typical velocity dispersion components, deduced by use of the virial 
theorem (hereafter referred to, in general, as the virial velocity dispersions) 
are global quantities, related to the virial potential energy of the subsystem as 
a whole, and so, by construction, independent of the specific orbital distribu- 
tion of the particles. This important property, however, is also a weakness of 
the virial theorem, in dealing with velocity dispersion components measured 
in the central region of a galaxy. In fact, it is well known that the related val- 
ues can be significantly different for structurally identical subsystems (and so 
characterized by identical virial velocity dispersion components), due to dif- 
ferent orbital structures (e.g., de Zeeuw & Franx 1991). When using central 
velocity dispersion components, an approach based on Jeans equations (even 
though still questionable) is to be preferred (e.g., Ciotti et al., 1996; Ciotti & 
Lanzoni 1997; Loewenstein and White, 1999). On the other hand, a compar- 
ison between the results obtained by use of either above mentioned methods, 
may provide additional support to both of them and/or useful indications on 
the nature of the problem under investigation. 

Strictly speaking, the central velocity dispersions (along the line of sight) 
in EGs, which are deduced from observations, should be scaled to the virial 
velocity dispersions. Both observational evidence (e.g., Gerhard et al. 2001) 
and theoretical arguments (e.g., Nipoti et al. 2002) point towards the exis- 
tence of dynamical homology in EGs. In particular, a linear relation is found 
between a local parameter, averaged central velocity dispersion, and a global 
parameter, inferred maximum circular velocity, cto.i = {'^/'i){vc)max (Gerhard 
et al. 2001). Accordingly, the central velocity dispersion components are ex- 
pected to be proportional to the virial velocity dispersion components. Then 
it could be assumed that the related proportionality factor is of the order of 
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unity. 

In fact, typical peculiar velocity component distributions within EGs 
show a maximum which is rapidly attained in the central region (at about 
1 kpc), and a slow decrease occurs moving outwards (no more than about 
13% the maximum at about 10 kpc), at least in the case of isotropic orbits; 
for further details see related attempts (e.g., Loewenstein and White, 1999). 
Accordingly, both the central and the virial velocity dispersion components 
are expected to be of comparable order, slightly less than the maximum of 
the peculiar velocity component distribution (e.g., Cappellari et al., 2006). 
On the other hand, most EGs arc moderately radially anisotropic (e.g., Ger- 
hard et al. 2001), and the related variation in central velocity dispersion (an 
increase for increasing af and vice versa) is also expected to be moderate. 

5.2 Input parameters, specific assumptions, and re- 
sults 

The main assumptions of the current model are (i) homeoidally striated 
density profiles and (ii) similar and similarly placed boundaries. Aspherical, 
hetcrogcncuos, self-gravitating fluids in dynamical or hydrostatic equilibrium 
exhibit isopycnic surfaces different from ellipsoids (e.g., Chandrasekhar, 1933; 
Chandrasekhar and Lebovitz, 1962; Vandcrvoort and Welty, 1981; Lai et al., 
1993), and the above assumption (i) is due to reasons of simplicity. 

Axisymmetric configurations with nonsimilar boundaries have been in- 
vestigated in the special case of homogeneous density profiles i.e. UU macro- 
gases, where the shape of one component is kept fixed and an additional 
variable, the axis ratio of the other component, together with an additional 
relation, the angular momentum conservation of the other component, must 
be considered (CS90). A main feature is that the isofractional mass curves 
on the {Oy(f)) plane cannot converge to at y — > -|-oo, as the ending point 
occurs when the inner subsystem attains a fiat configuration. For further 
details refer to the parent paper (CS90). Being the calculations for hetero- 
geneous density profiles much more complicated, the above assumption (ii) 
is also due to reasons of simplicity. It can be expected that the effect due 
to nonsimilar boundaries is maximum for homogeneous density profiles, and 
decreases as the profile is steeper, to be null for mass points surrounded by 
a massless atmosphere (Roche ellipsoids). Then the results found for homo- 
geneous density profiles (CS90) make a valid reference for inhomogeneous 
density profiles. In conclusion, the above assumptions (i) and (ii) are related 
to a viable model which can be used for specific applications to large-scale 
celestial objects, in particular EGs. 
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Given a typical EG, a natural question (in the light of the model under 
discussion) arises about its position on the (Oycj)) plane for assigned density 
profiles and specified input parameters. Towards this aim, the following 
main assumptions are made: (a) the stellar and the DM distributions are 
described by homeoidally striated, similar and similarly placed, HH or HN 
density profiles, respectively; (b) undetected baryons trace either DM haloes 
or EGs; and (c) the virial theorem holds for each subsystem within the tidal 
potential of the other one. 

According to the above assumptions, a typical EG embedded in a DM halo 
is idealized as two homeoidally striated, similar and similarly placed matter 
distributions, where the star and non baryonic subsystem are described by 
HH or HN density profiles, respectively. 

For assigned density profiles, the macrogases equation of state, Eq. fpUl) . 
depends on two independent variables: the fractional truncation radius, y = 
Rj/Ri, and the fractional mass, m = Mj/Mi, or their scaling counterparts, 

= rj/rj and = M^/M}. To represent celestial objects in the {Oycp) 
plane, two additional equations are requested. Having in mind an application 
to EGs and their hosting haloes, denoted in the following by the index i 
and j, respectively, for y > 1, an additional relation between the stellar 
projected velocity dispersion, {(Ti)R^/s, averaged over the aperture used for 
spectroscopic observations (r = Re/ 8), and the fractional scaling radius, 
y\ may be deduced from the virial theorem related to the star subsystem, 
Eq. ( I8b1) . by solving the appropriate Jeans equations, following an earlier 
approach (Ciotti et al., 1996). The result may be written as: 

(^i)L/8 = ^T^V'^(l/^"^^) ; (75) 

where ipi is a function calculated by a numerical algorithm. 

Using Eq. (ffSll implies further restrictions, as its validity is limited to (i) 
spherical-symmetric matter distributions; (ii) isotropic peculiar velocity dis- 
tributions; (iii) infinite truncation radii; (iv) HH density profiles. Concerning 
points (i)-(iii) mentioned above, acceptable approximations may safely be ex- 
pected. On the other hand, HH density profiles provide a viable description 
to EGs embedded within DM haloes (e.g., Ciotti, 1996). 

At this stage, an additional relation is needed. The mere existence of 
a fundamental plane (Djorgovski and Davis, 1987; Dressier et al., 1987) in- 
dicates that structural properties in EGs span a narrow range, suggesting 
that some self-regulating mechanism must be at work during formation and 
evolution. In particular, projected light profiles from EGs exhibit large de- 
gree of homogeneity and may well be fitted by the r^/"^ de Vaucouleurs law. 
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Accordingly, a narrow range may safely be expected also for EG fractional 
mass and the assumption, m = const, appears to be a viable approximation. 
The last, together with Eqs. (l26ll and ( TfSj) . makes a system of three equations 
in the three unknowns, 0, y, m, via Eqs. (1211) . for HH density profiles. Then 
the position of EGs and their hosting haloes on a selected isofractional mass 
curve in the {Oyip) plane, can be determined. 

The sample used {N = 16) is extracted from a larger sample (A^ = 25) of 
EGs and lenticular galaxies investigated within the SAURON project (Cap- 
pellari et al., 2006, Tab. 1 therein), for which the parameters of interest i.e. 
masses, projected central velocity dispersions, and effective radii, can be de- 
rived from the data listed in Tab.[2j More specifically, with regard to the 
inner (stellar) subsystem, projected central velocity dispersions are deduced 
from the luminosity-weighted second moment of the line-of-sight velocity 
distribution within the effective radius, cxe, as (cTi)ij^/8 = 8°°^^(Te; masses 
are deduced from luminosities and mass-luminosity ratios (in J-band), as 
M,/Mio = (L/L0)[(M,/L)/(lOiOm0/L0)]; L/Lq = exp,,{-OA[lT - (m - 
M) — 4.11]}; scaling radii are calculated from effective radii (in arcsec) and 
distances, by use of a profile factor, equal to 1.81, related to the case un- 
der discussion (Hernquist, 1990), as r|/kpc = (_Re/kpc)/1.81; Re/kpc = 
[(i?e/arcsec)(d/Mpc)]/206.265; d/Mpc = expio[(m - M)/5 - 5]. For further 
details refer to the parent paper (Cappellari et al., 2006). 

The substitution of the selected value of the fractional mass, m, into 
Eq. (1751) . allows the value of the fractional scaling radius, y^, for each sample 
object. Finally, the substitution of m) values into Eq. ( !26l) . particularized 
to HH density profiles via Eqs. ( I56l) -( l59|) . allows the value of the fractional 
virial potential energy, 0, and the position of each sample object on a selected 
isofractional mass curve in the (O?/0) plane, can be determined. 

The results are plotted for different choices of scaled truncation radii, 
(Sj, Hj), and different choices of fractional masses, m, in Fig. HI Isofractional 
mass curves, m = 10 (top panels), and m = 20 (bottom panels), are labelled 
by the selected choices of scaled truncation radii, (Sj, Sj), indicated near the 
corresponding curves, and the sample objects are represented as dots. Curves 
lie above and below the critical curve for m = 10, while all curves lie above 
the critical curve for m = 20, with respect to the cases considered. 

The above description is due to reasons of simplicity, but the model does 
not necessarily imply that sample objects must be located on the same 
isofractional mass curve and/or correspond to the same scaled truncation 
radii. What is relevant is the position of EGs on the (O?/0) plane, which 
implies the following assumptions: (i) there are in the universe 16 EGs with 
intrinsic values of the parameters equal to their counterparts listed in Tab. [2] 
for sample objects, and (ii) related star distributions are well described by H 
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Table 2: Data related to a subsample (A^ = 16) extracted from a sample 
(A^ = 25) of EGs and lenticular galaxies investigated within the SAURON 
project (Cappellari et al., 2006), and deduced values of the parameters of 
interest. Column caption: (1) NGC number; (2) effective (half-Hght) ra- 
dius, Re, measured in the 7-band; (3) total observed /-band galaxy magni- 
tude; (4) mass-luminosity ratio (including DM) deduced from the best fitting 
three-integral Schwartzschild model, computed at a fiducial inclination; (5) 
mass-luminosity ratio of the stellar population; (6) galaxy distance modu- 
lus (hats avoid confusion with the fractional mass, m, and the total mass, 
M); (7) galaxy mass, calculated as Mi/Mio = (L/L0)[(Mj/L)/(lO^°m0/LQ)]; 
L/Lq = expiQ{—OA[lT — {rh ~ M) — 4.11]}; (8) galaxy central velocity dis- 
persion, deduced from the luminosity-weighted second moment of the line- 
of-sight velocity distribution within the effective radius, a^, as {cFiju^/^ = 
g0.066^^. ^g-j galaxy scaling radius, calculated as r|/kpc = (i?e/kpc)/1.81; 
i?e/kpc = [(i?e/arcsec)((i/Mpc)]/206.265; d/Mpc = expio[(m - M)/5 - 5]. 
The factor, 1.81, is related to an assumed Hernquist profile for the inner 
subsystem (Hernquist, 1990). The factor, 206.265, is related to the choice of 
measure units. For further details refer to the parent paper (Cappellari et 
ah, 2006). 
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density profiles. 

In fact, observational uncertainties on the quantities of interest (Cappel- 
lari et al., 2006) make large errors on the fractional scaling radius, y'^ . As an 
example, a single galaxy, NGC 3379, and a single parameter, the effective ra- 
dius, shall be considered. In the case under discussion. Re = (42.0 =F 7.1) arc- 
sec (Cappellari et al., 2006), but different estimates exist, such as (54.8 =F 3.5) 
arcsec (Capaccioli et al., 1990) which is consistent with the above result 
within 2aR^. Using the latter value, the fractional scaled radius, y^, passes 
from 5.74 to 4.08 for a fractional mass, m = 10, and from 8.45 to 6.07 for 
m = 20. 

Different values of parameters, m, Sj, Sj, makes sample objects locate on 
different isofractional mass curves, as shown in Fig. HI where y depends on 
m but is independent of Sj and S^. Fiducial values of the above mentioned 
parameters, say 10 < m < 20, 5 < H„ < 20, u = i,j, define a region on 
the {Oycp) plane, approximately as: 2.5 < y'^ < 6.0, 7 < < 11, m = 10; 
4 < 1/^ < 9, 17 < < 23, m = 20; where sample objects are positioned, 
as shown in Fig.O Dots and asterisks correspond to inner scaled trunca- 
tion radii, Sj = 5 and 10, respectively. Small and large symbols represent 
sample objects and configurations where the virial potential energy of the 
inner subsystem attains the maximum value with respect to a frozen outer 
subsystem, respectively. No such configuration exists in the cases considered, 
for fractional masses, m ~ 20. The related parameter space is restricted to 
a region close to the angle defined by the locus of large dots and asterisks, 
respectively. The intersection between the above mentioned loci occurs for 
outer scaled truncation radii, Sj, slightly larger than 10. The change of po- 
sition for NCG 3379, due to a change in effective radius, from Re = 42.0 to 
Re = 54.8, is indicated by squares, regardless of the vertical scale. 

In the special situation where the outer subsystem remains frozen and 
the inner one is free to contract or to expand, with unchanged mass and 
density profile, the virial potential energy of the inner subsystem may attain 
an extremum point of maximum (Secco, 2000, 2001, 2005; Marmo and Secco, 
2003). For HH density profiles, lower fractional masses (m = 10) yield no ex- 
tremum point, while the contrary holds for larger fractional masses (m = 20). 
The special configurations related to the maximum virial potential energy of 
the inner subsystem, when the outer subsystem remains frozen (hereafter 
quoted as "the maximum configuration"), are represented as large dots in 
Fig.m (bottom panels) and large dots and asterisks in Fig. [51 

The maximum configuration appears to have little relevance in the light 
of the current model for a number of reasons. First, it occurs for fractional 
masses above a threshold. Second, the loci of maximum configurations on 
the {Oycf)) plane show no correlation with the parameter space of the model 
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related to sample objects. More precisely, the loci of maximum configurations 
are narrower and extend from bottom left to top right, while sample objects 
lie on a broader and less inclined band, as shown in Fig.O 

The position of the maximum configuration on the isofractional mass 
curves, m = 20, depends on the concentration of the outer subsystem, while 
the position of sample objects can only be vertically shifted, as the fractional 
scaling radius, |/^, is independent of the concentration, via Eq. (173|) . In this 
view, it is difficult to conceive maximum configurations as ending points of 
any evolutionary track on the {Oycj)) plane, even if observational uncertainties 
are high. 

EG position on the {Oycj)) plane, in particular along an isofractional mass 
curve (m = const), represents the ending point of evolutionary tracks on the 
above mentioned plane. The related configurations may be thought of as 
virialized to a good extent, in that sample objects listed in Tab. [2] show no 
sign of ungoing merger and star formation burst. The coincidence of ending 
points with maximum configurations would imply, for fixed (m, Sj,Sj): (i) 
a single ending point for all sample objects i.e. homologous evolution, and 
(ii) m above the threshold for the occurrence of maximum configurations, 
which does not necessarily happen, as shown in Fig. HI On the contrary, the 
evolution of EGs appears to be non homologous, in the sense that different 
ending points take place along the selected isofractional mass curves defined 
by the HH macrogases equation of state. 

Owing to Eqs. ^ and ffT2bl) . the total energy of the virialized system 
reads: ^ 

^ = + (-^j*)vir] = — (-E'j)kin " (-^jOkin ; (76) 

and the combination of Eqs. fll2b|) . fl2^ . and flTUj) yields: 

E = \{1 + <Pm,U = -(1 + <l>mhn ; = ^fl^ ; (77) 

where the kinetic energy of the inner subsystem, (i?j)kin, may be deduced from 
observations, and the related total energy, E, may be read on the vertical 
axis of Fig.m Within the range, (pram < < 4>max, three configurations exist 
with same virial potential (or kinetic) energy ratio, 0, and fractional mass, 
m, but different fractional truncation radius, y. 

The restriction to constant fractional mass also in time, makes evolution- 
ary tracks on the {Oycj)) plane locate on the related isofractional mass curve, 
(m = const), where displacements from the left to the right (increasing y val- 
ues) are due to energy dissipation, and displacements from the right to the 
left (decreasing y values) are due to energy acquisition. In the former alter- 
native, changes in and (i?jj)vir must act to yield a decreasing E (increasing 
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in absolute value), the larger energy change being related to the larger y 
change. For further details refer to Appendix [Bl In this view, the maxi- 
mum configuration corresponds to a special energy change, (Ai?)niax, and to 
a special y change, (A?/)max- But there is no apparent reason for which a 
special amount of energy has to be dissipated starting from the beginning 
of evolution, to attain the maximum configuration. A similar result holds in 
the latter alternative. In fact, it can be seen in Fig.lHthat the sample objects 
show no connection with the maximum configuration. Then the maximum 
configuration, which implies a frozen outer subsystem, has no special rele- 
vance in the light of the current model. For an analysis of different theories 
on the maximum configuration, refer to a specific study (Valentinuzzi 2006, 
Chaps. 3-5). 

The above results hold if the baryonic subsystem is mainly in form of stars. 
Let us take into consideration a different scenario, where a less concentrated 
gaseous subsystem than the stellar one is also present, as in the Coma cluster 
of galaxies (e.g., Briel et al., 1992), and assume the same mass distribution 
as in the non baryonic matter to preserve use of two-component models 
(Caimmi, 2003). Accordingly, a typical EG is idealized as formed by an 
inner subsystem made of stars and an outer subsystem made of gas and non 
baryonic matter. As the amount of baryonic and non baryonic matter have 
to remain unchanged, the inner and the outer subsystem are less and more 
massive, respectively, than in absence of undetected baryons. Again, it is 
assumed that the related mass distributions are represented by HH density 
profiles. Then a similar procedure may be followed, keeping in mind that the 
fractional mass, m, is the ratio of nonbaryonic + extragalactic gas mass to 
star mass. For further details refer to an earlier attempt (CM03). 

6 Conclusion 

Two-component systems have been conceived as (two-component) macro- 
gases, and the related equation of state has been formulated using the virial 
theorem for subsystems (Limber, 1959; Brosche et al., 1983; Caimmi et al., 
1984; Caimmi and Secco, 1992), under the restriction of (i) homeoidally stri- 
ated ellipsoids (Roberts, 1962) and (ii) similar and similarly placed bound- 
aries. 

Explicit calculations have been performed for a useful reference case and 
a few cases of astrophysical interest, both in presence and in absence of 
truncation radius. More specifically, the following cases have been dealt with: 
IJ= UU, PP, HH, HP, HN, where I and J denote the inner and the outer 
density profile, respectively, and the other captions relate to the following 
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density profiles: U (p =const), P (Plummer, 1911), H (Hernquist, 1990), N 
(Navarro et al., 1995, 1996, 1997). Shallower density profiles (UU, PP), have 
been found to yield an equation of state, = (f){y,m), characterized by the 
occurrence of two extremum points, one maximum and one minimum, as 
in an earlier attempt (CS90). Steeper density profiles (HH, HP, HN), have 
been found to produce a similar equation of state where, in addition, a single 
horizontal inflexion point occurs in a critical isofractional mass curve, and 
isofractional mass curves related to lower values, m = Mj/Mi < merit, show 
no extremum point. The similarity between isofractional mass curves and 
van der Waals' isothermal curves, has suggested the possibility that a phase 
transition could take place in a bell-shaped region of the {Oycf)) plane. 

Further investigation has been devoted to HH density profiles for which 
a numerical algorithm (Ciotti et al., 1996) has been used to represent EGs 
and their hosting DM haloes along selected isofractional mass curves on the 
{Oy(j)) plane, under the assumption that the related fractional mass has the 
same value in different systems. In the light of the model, the evolution of 
EGs has been found to be non strictly homologous, in the sense that the end 
of evolutionary tracks on the {Oyip) plane occur at different points along the 
related isofractional mass curve, instead of being close to a single point. 
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Appendix 

A Plummer density profiles and n = 5 poly- 
tropes 

The Lane-Emden equation reads (e.g.. Lane, 1870; Emden, 1907; Chan- 
drasekhar, 1933; Horedt, 2004, Chap. 2, §2.1): 

' ' '&.^] - -r ; (78a) 



ClE d^LE V~'^^ '^'^LE , 

^(0) = 1 ; ^(Hle) = ; (78b) 
r R 

^LE = — ; Sle = — ; (78c) 

OLE «LE 

p(r) = Ar(eLE) ; (78d) 

where n is the polytropic index (0 < n < 5 and 0.5 < n < 5 for realistic 
collisional and collisionless equilibrium configurations, respectively), ole is a 
scaling radius and A is the central density. 

In the special case, n = 5, the integration of Eq. fl78al) yields (Schuster, 
1883; see also Chandrasekhar, 1939, Chap. IV, §4; Horedt, 2004, Chap. 2, 
§2.3.3): 

/ 1 \-i/2 

^(^le) =(l + -^leJ ; Sle -> +00 ; n = 5 ; (79a) 
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.o d^ 



-3/2 



LE 



.73, 



'3/2 



lim 

6le-^+oo 



"^LE' 



d^ 



LE, 



(79b) 
(79c) 
(79d) 



and the general expression of a polytrope mass (e.g., Chandrasekhar, 1939, 
Chap. IV, §4; Horedt, 2004, Chap. 2, §2.6.3): 



d9 



d^: 



(80) 



LE / 



in the case of interest, using Eqs. (]79l) . reduces to: 

M = VSATiXal^ ; n = 5 ; (81) 

finally, the explicit expression of the density profile results from the combi- 
nation of Eqs. ( 178dl) and f l79al) . as: 



p(r) = A ( 1 + -^LE 



1 \ -5/2 



n = 5 



(82) 



and the following relation is derived from comparison with Eqs. fll4al) and 
(ETD: 



which, in turn, implies the following: 

^LE ^ ■ 



A ( 1 + g^LE 



1 \ -5/2 
'2 



7=0; a =2 
A = 25/V ; 



5 /?-7 /5 



2 ' 



and the substitution of Eqs. fl84ap and fl84cl) into (IH^ yields: 

pW = pV(0 ; 

95/2 



(1 + ^2)5/2 



(83) 

(84a) 

/3 = 5 ; (84b) 
(84c) 

(85a) 
(85b) 



according to Eq. (1371) . 
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An equivalent formulation can be obtained by the combination of Eqs. (1811) 
and (IH^ . The result is: 

, , 3M (rt)2 
which is known as the Plummer (1911) density profile. 



B Quasi static contraction in presence of tidal 
potential 

The scalar virial equations, Eqs. ( 112bp . for assigned density profiles i.e. fixed 
f^ui 7«, and S„, u = depend on four parameters e.g., masses and 
scaling radii, and r^. The further assumption of spherical-symmetric 
matter distributions is only to simplify calculations. Let an assigned amount 
of energy, AE < 0, be instantaneously dissipated within the inner subsystem, 
as: 

(Ei)kin ^ (^i)kin + ; (87) 

and the system readjust to attain a new configuration in accordance with 
Eqs. (112bp where, in general, density profiles and scaling radii are changed, 
while masses and the outer subsystem truncation radius may be assumed 
fixed, or their variations assigned. More specifically, density profiles are 
steepened by energy dissipation, and the system gets bounder i.e. a larger 
(in absolute value) binding energy is attained. 

The further constraint of density profiles unaffected by energy dissipation, 
implies the following changes {u = i,j): 

^M, + AM„ ; AM, = ; (88a) 

S„-^H„ + AS„ ; AS„ = ; (88b) 

Ru^Ru + ARu ; ARj = ; (88c) 

rl^ri + Ari ; Ar] = ; (88d) 

(^«)kin ^ (^«)kin + AE, ; AEi + AEj = AE ; (88e) 



where the assumption that energy dissipation within the inner subsystem has 
no effect on (i) masses, M,; (ii) scaled truncation radii, S,; density profiles, 
{au, Pu,1u)', (iv) truncation radius of the outer subsystem, Rj; implies a 
frozen outer subsystem. In general, the changes, AM, = Cm^u, Ai?, = 
CuRu, and Arl = (u^l, which make m and conserved, should be specified. 
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The application of the scalar virial theorem, Eqs. (1121) . to the subsystems 
before and after energy dissipation, reads: 

?7,,(yt) + 2(Ejkin = ; (89a) 
[/„,(l/^ + A|/t) + 2[(E„)kin + AE„] = ; (89b) 
Uuv = {Euv)vir ; u = i,j ; v = j,i ] (89c) 

where the fractional scaling radius, = rj/rj, has been chosen as variable. 
The combination of Eqs. (I89ap and (189bl) yields: 

A{Euv)vir + 2AEu = ; u = i,j ; v = j,i ; (90a) 
A(^™)vir = t/„,(i/^ + A|/t)_t/„,(i/^) ; (90b) 

and the combination of Eqs. f l90al) . by use of f l88ep produces: 

A{Ei,U, + A{Eji),i, + 2AE = ; (91) 

which is a transcendental equation in A?/^ provided the density profiles and 
the amount of dissipated energy, AE, are specified. In general, the changes, 
AMu = CmMu, ARu = CuRu, and Ar^ = Cu^l, which make m and Eu con- 
served, should also be specified. Then the remaining parameters related 
to the relaxed system, due to energy dissipation, may be determined. In 
particular, Ay"^ > is expected together with m = const due to mass 
conservation in each subsystem or mass variation of the kind considered, 
m = {Mj + AMj)/{Mi + AMi) = Mj/Mi. Accordingly, the position of the 
system in the {Oycf)) plane moves from the left to the right along the related 
isofractional mass curve (see e.g., Fig.H]). 
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Figure 1: Isofractional mass curves related to (from top left in clockwise 
sense) UU, PP, HH, HP, macrogases for S —>■ +00, 'E^j/Ei 1. Thick curves 
represent the locus of minimum (left branch) and maximum (right branch) 
points. The critical curve is also thickened, and the horizontal inflexion point, 
or critical point, is defined (when present) as the tangent point with respect 
to the locus of the extremum points. The value of the fractional mass is 
m = 0(1, m) and can be read on the intersection between the selected curve 
and the dashed vertical line, y'^ = y = 1, provided the density profile of the 
inner and the outer subsystem belong to the same family (UU, PP, HH). 
All the curves diverge at y ^ ancl^converge to at y — > +00. For UU 
macrogases, the extremum points of the m = isofractional mass curve occur 
at y = and y = respectively, and the locus of the minimum points 

has a vertical asymptote, y = \/5/3. For PP macrogases, the extremum 
points of the m = isofractional mass curve occur at y = and y ^ 1.20, 
respectively, and the locus of the minimum points has a vertical asymptote, 
y ^ 0.84. 
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Figure 2: Isofr actional mass curves related to HH macrogases, for different 
choices of scaled truncation radii, (Sj, Sj), represented on each panel, where 
(from bottom to top) m =10, 20, 60. The critical isofractional mass 
curve, where the horizontal inflexion point, or critical point, is marked by 
a St. Andrew's cross, is thickened on each panel. All the curves diverge at 
y — > and converge to at |/ — > +oo. 
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Figure 3: Isofr actional mass curves related to HN macrogases, for different 
choices of scaled truncation radii, (Sj, Ej), represented on each panel, where 
(from bottom to top) m =10, 20, 60. The critical isofractional mass 
curve, where the horizontal inflexion point, or critical point, is marked by 
a St. Andrew's cross, is thickened on each panel. All the curves diverge at 
y — > and converge to at y — > +oo. 
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Figure 4: Isofractional mass curves, m = 10 (top panels), and m = 20 (bot- 
tom panels), related to HH macrogases, for different choices of scaled trun- 
cation radii, (Sj,Sj), indicated near the corresponding curves, and related 
positions of = 16 EGs (dots), listed in Tab. El Different cases correspond 
to vertical shifts of the sample objects. Larger dots represent configurations 
where the virial potential energy of the inner subsystem attains the maximum 
value with respect to a frozen outer subsystem. 
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Figure 5: The parameter space of the cases plotted in Fig. HI roughly out- 
lined by the shaded region. Dots and asterisks correspond to inner scaled 
truncation radii, Sj = 5 and 10, respectively. For fixed fractional mass, m, 
the corresponding outer scaled truncation radii read = 5, 10, 20, from 
top to bottom. Small and large symbols represent sample objects and con- 
figurations where the virial potential energy of the inner subsystem attains 
the maximum value with respect to a frozen outer subsystem, respectively. 
No such configuration exists in the cases considered, for fractional masses, 
m ~ 20. The related parameter space is restricted to a region close to the 
angle defined by the locus of large dot^^and asterisks, respectively. The inter- 
section between the above mentioned loci occurs for outer scaled truncation 
radii, Ej, slightly larger than 10. The change of position for NCG 3379, due 
to a change in effective radius, from Rg = 42.0 to Re = 54.8, is indicated by 
squares, regardless of the vertical scale. 



